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O ' Abstract 

o 

We consider departures from hamiltonian dynamics in the evolution of neutral kaons 
due to their interactions with environment that generate entanglement among them. 
We propose a phenomenological model of stochastic re-scattering and estimate the 
coefficients of the effective hamiltonian and decoherence terms. Finally we analyze the 



in 



interplay between the weak interaction and matter effects and propose the observables 



£Sj | suitable to measure the matter characteristics of CPT violation and decoherence. 

> 

<N 1. Introduction 

o ■ 

O Since the fifties the system of neutral kaons has provided us with a unique laboratory to 
investigate the superposition principle, interference effects and CP violation |[|, |2[. CP 
violation was first measured in || long ago, in the mixing of the neutral kaons, the so called 
^P" 1, indirect CP violation, whereas it was not until recently that direct violation, in the decay 
amplitudes to two pions, was experimentally established (see 0, |J, ||). There are at 
present a few dedicated experiments in project to make precision measurements of the CP 
>- \ breaking parameters ||, ||, [0. The analysis of their data may require to take properly 
into account the effects of decoherence due to entanglement of the neutral kaon system with 
the environment in which it evolves, which entails a pure kaon state to convert into a mixed 
one. We consider the evolution of neutral kaons in the presence of matter, in an environment 
that is not the perfect vacuum, which generates departures from hamiltonian dynamics in 
the kaon sub-system. These are effects that exist in addition to the weak interactions and 
are dominated by the strong interactions of the kaons with the environment: they lead to an 
effective breaking of charge conjugation, because the environment is made of matter, not of 
anti-matter; and to an effective breaking of time-reversal invariance associated to the large 
number of degrees of freedom in the environment that leads to irreversibility. 

Similar analysis can be found in the literature ||, |§, p0| , ||11|| . There, the motivation is 
to study decoherence effects that come from quantum gravity. Ours is the interaction with 
the environment that is not the vacuum HT2I. 
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In section 2 we explain how the propagation in an environment being projected to the 



kaon system is described by a Lindblad master equation |13| which characterizes the de- 
parture from hamiltonian evolution; we introduce a criterion to quantify decoherence and 
build a phenomenological model based on a stochastic re-scattering of kaons in a dense 
matter induced by strong interactions. In section 3 we restrict the general analysis to a 
two-component K° — K° system and estimate coefficients of both effective hamiltonian and 
Lindblad interaction. In particular we conjecture that the decoherent part may be equally 
important for both hamiltonian and Lindblad piece. In section 4 we analyze the interplay 
between the weak interaction and matter effects and propose the observables suitable to 
disentangle the parameter responsible for CPT violation in matter 5k, the characteristic of 
indirect CP violation in vacuum e, and the decoherence parameter X33. In the appendix we 
propose a toy model of two level environment to explain how the Lindblad interaction comes 
about, in the Markovian approximation. 

The neutral kaon evolution in the vacuum is analyzed with the Schrodinger equation 

ip{t) = H w p{t) - P {t)Hl (1) 



where H w is the LOY [14] effective hamiltonian H w = ^2 I=S L XAK)(Ki\, which is not 
hermitean and its eigenstates Xj = mi — -7/ possess a non-vanishing imaginary part of 
their decay rates. S, L stand for the short and long components, and \Kj) are such that 
(Ki\Kj) = 5 U . Recall pf that m s ~ m L ~ 500 MeV, Am = m L -m s = 3.5 x 10~ 12 MeV, 



75 = 7.3 x 10~ 12 MeV and 7 L = 1.3 x 10~ 14 MeV. 

The \Kj) have diagonal evolution \Kj(t)) = e~ tXjt \Kj(0)) . In terms of the CP eigenstates 

\Ks) = 1 1 i\Ki) + e s \K 2 )) , \K L ) = 1 (\K 2 ) + e L \K x )) ; (2) 

+ M 2 \A + M 2 



they reduce to the CP eigenstates \K±) = \ (\K°) + \K )), \K 2 ) = \ (\K°) - \K )), if CP 
is conserved. 

At this point we recall that CPT conservation requires = e^; T invariance, eg + e^ = 0; 
and CP, e s = e L = 0. 

Indirect CP violation has been measured |§, fl5[ 



e ~ 2.3 x 10- 3 e in/4 , (3) 
(CPT conservation is assumed, = ex = e.) Direct CP violation has also been measured 

e'/e ~2x 10" 3 , (4) 

where 

A(K L -> 7T+7T-) _ , A(K L -> 7r 7T ^ 



e L + e j "7777 rTTTT = e L ~ ^ ,^L 



A(K S -> 7T+7T-) u ' A(K S -> 7r°7r°) 



e. 



2. Propagation in environment 

The propagation in an environment is dominated by the strong interactions, and we postpone 
to the next section the inclusion of weak interactions. The evolution of the total system 
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[Kaon + (Large) Environment] is unitary, with the time evolution operator given by U(t) = 
exp(-iHxotait)', whereas the dynamics of the kaon sub-system alone is obtained by tracing 
out the environment degrees of freedom: 

p K {t) = Tr E NV (U(t)p 

K+ENV 

(O)tft(t)) . (5) 

This gives a complicated evolution for px, even if the initial conditions are a direct product 
Pk+env(Q) — Pk(®) ®Pbw(0)- However, when the interaction is tiny (e.g., as in the case of 
a diluted environment) the dynamics of px becomes approximately free from memory effects 
[|T6f (Markovian), and it is dictated by a Lindblad master equation |TJj of the form 



p = -i (H eff p(t) - p{t)H\ }f ) + L[p], (6) 

where 

L[P) = E ( A M] - \pA\A, - \a)A oP ) , (7) 

with the effective hamiltonian being in general non-hermitian. (Henceforth p stands for 
Pk)- Memory effects are being neglected, (see a toy model in the appendix). Notice that 
the new piece is linear in p but quadratic in the unspecified operators Af the Lindblad 
equation encodes a new dynamics that is not purely hamiltonian. The new term may entail 
decoherence by making pure states evolve into mixed states. Since the opposite process does 
not occur, i.e, no mixed state evolves to a pure one, we demand that the von Neumann 
entropy — plogp should not increase, under evolution with the new piece only. 

At this point we introduce a criterion to quantify decoherence. In its time evolution, 
there are two distinct effects that will affect the density matrix, namely the loss of quantum 
coherence, due to entanglement with the medium, and probability depletion, due to decay of 
the initial states which decreases tr (p). It is natural to introduce a new variable p through 
the relation p = p tr (p) , so that tr (p) = 1 remains normalized to unity. Notice that this 
conditional p matrix verifies the following closed evolution equation 

p = -i(H eff p - pH J eff ) + L[p] - 2tr (p lmH eff ) p. 

In terms of p we propose the quantity d = 2 detp = (1 — tr (p 2 )) > as an order parameter 
for decoherence: only pure states have d = 0. This parameter represents a linear entropy for 
the conditional p matrix [17j . 



Let us consider the particle scattering in a dense matter and derive an average description 
of its propagation for a relatively large, macroscopic time of flight. First we introduce the 
elementary evolution of the density matrix p induced by the scattering matrices Sj and as- 
sume a stochastic property of the medium, effectively described by the classical probabilities 



p{r + Ar) = Y.r ] S J p{r)S). (8) 

j 

We suppose that all possible scattering processes are included and the description is complete 
so that 

5>i = l. (9) 
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Thus for the unitary (elastic) evolution, £• 1 = Sj one obtains the conservation of probability: 

tr(p(r + Ar))=tr(p(r)) =1. (10) 

Further on a certain inelasticity of scattering can be taken into account by S^ 1 ^ Sj and 
\\SjS}\\ < 1. 

In order to implement the smooth average evolution the time intervals At should be 
chosen of macroscopic size I which is much larger than the size of a cristal link or a distance 
between molecules in a gas or a liquid. Therefore At ~ l/v ^> Rmt/v for a given particle 
velocity v and a radius of strong interaction Ri n t- 

Let us introduce the scattering amplitudes Sj — I + iTj. Then the one-step evolution 
can be characterized by Ap = p(r + At) — p(r) in a differential form: 

^ = -^^p-pHl^+^^^-liT^p + pT^)) (11) 

which represents a Lindblad equation for a coarse-grained time evolution. Herein the effective 
hamiltonian is given by 

Heff = J2^-(-T 3+ t -Wi). (12) 

This definition leads to the hermitian hamiltonian in the elastic case (TjTj = 2ImTj, a 
relation from unitarity) 

H eff = -J2%-KeT j = Ht ff . (13) 

3 

Respectively in (^|) the Lindblad matrices are Aj = \f^Tj. 

We see that the unitary evolution in a time-slicing approximation is generated not only 
by a hermitian hamiltonian H e ff but also by the unavoidable Lindblad interaction which 
does not necessarily entail the decoherence effects. The latter ones arise only if the medium 
is stochastic with Tj ^ 1. 



3. Kaons in a medium 



The operators Aj are governed by the strong interactions and the properties of the medium. 
As to the hamiltonian we have to sum up the contributions that come from the original 
evolution in the vacuum (weak interaction ([!])) plus the rescattering effects already contained 
in (|i"2"|) H = H w + H e ff-. therefore the constants ex,, es are called el, e~s and contain the effects 
due to evolution in a medium of matter. Since strong interactions conserve strangeness, which 
in the absence of the weak interactions would become a superselection quantum number, we 
further demand that any density matrix, at large times, should end up as a mixture of \K°) 
and \K°). The most general parametrization that satisfies the above considerations leads, 
for the two state neutral kaon system, to a non-hamiltonian contribution in (Q) 



Hp) 



(° 








^ 




f p° 


\ 





Si 










P 1 










S22 


$23 




P 2 




V o 





S23 


S33 J 




VP 3 





(14) 
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where in the CP eigenbasis p = \{p°I + p- a) (see (0)). We have included a piece in H e ff 
from L[p] which is of hamiltonian form and left in L[p] the genuine non-hamiltonian part. 

The matrix is symmetric, Si, S22, S33 > 0. They also verify the relations ||11|| , providing 
the complete positivity, 

s l < $22 + S33, S33 > (Si 



°33 

S22 > ( s l 



S 2 2 



S33 



S22 < Si + S33, 
s 33 < s l + S22, Si > (S22 — ^33)^. 

Now let us find the structure of the evolution equations for the two-level system of kaons 
(K ,K ). The conservation of strangeness implies that the S-matrices take the following 
form 











where Sj 



(±) 



2 i s j -'- s j/ 



Respectively, 



T; 



4 w 



(15) 



(16) 



Accordingly the effective hamiltonian and the Lindblad interaction look as follows: 



H eff : 

L[p] = t 

we used the following notations, 

t± = 



-t + + - 2 (t 



+ t__) +<T 3 (-*- + **+-, 

(a 3 pa 3 -p), 



(±). 



(17) 



Thus, in our approach s 2 2 = s 33 = —It , s\ = s 2 3 = 0. There are other ||, [18| models 

where these cofficients do not obey this constraints, when the effects are induced by quantum 
gravity. 

Let's adopt the model of two-component evolution with r% — 1 — r-i, r 2 = O(At) <C 1. 
Respectively we suppose that ^1 -C fol- Then the coherent part is dominantly described 
by ti and its nontrivial contribution to the non-hamiltonian evolution ([17]) can be estimated 
as follows: 



= Ar 1 



tr (a 3 H ef f) I Ar < (Am) 2 10" 2 Ar, 



(19) 



where we have used the Good-Kabir theory for the coherent regeneration [19] in which the 
forward scattering is dominating: 



2nu 
m K 



f 



(±) 



with the forward scattering amplitudes /j ± ' ) ; 
per unit volume and being a kaon mass. 



\t±\ < Am ■ KT 1 (20) 
f K ± f K , v standing for the particle density 
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One can further estimate At taking into account that the optical picture arises at macro- 
scopic distances, say I ~ 10~ 4 cm as a collective effect of many scatterers. Then At — l/v~l, 
for relativistic kaons (v ~ 1) which entails AmAr ~ 10~ 5 and therefore L\ ~ Am 10~ 7 . This 
estimate can be considered as an upper bound for the deviation from a hamiltonian evolution 
due to averaging in time. It proves it to be difficult to observe these coherent effects. 

On the other hand, the second, decoherent part may equally contribute both in the 
effective hamiltonian and in the Lindblad coefficients. Indeed, if r 2 = O(At) <C 1 but 
|*i I *C 1^1 (impurities of heavy nuclei or density clusters) one may expect that \ti\ ~ r 2 |4 '| 
and further on 

At 1 At 1 1 

Thus the decoherent part of the Lindblad may be well observable. We nevertheless take here 
a conservative estimation for the Lindblad coefficients to be of order Ame|. 

4. Decay probabilities 

The complete hamiltonian of ( |i~7"D contains two pieces: the first, proportional to the iden- 
tity, affects the mass and the decay rates by a shift of a same amount; the second piece, 
proportional to 03, modifies the mixing parameters to ei = e mc + 8l, e~s = tvac — $l, leaving 
the masses and decay rates unaffected to leading order in the perturbation. Here 8i is the 
parameter that measures the CPT breaking due to matter effects, which is O(10~ 2 ) in a 
dense matter |B|. In principle, 5 can be measured (see 



The decoherence effects encoded in the Lindblad interaction modify the evolution pa- 
rameters as follows 

lL = lL + — , IS = IS + "Tp 7 = 75 + 1L + S 2 2 = 1L + IS 



s 33 S22 
£33 = , X = — - = /X33, 

2-fs 2 Am 



1 

A\ = \ L - As = (m L - ms) + -{is ~ 7l), 

where both mi — ms, 7s — 7l > 0. Herein the masses m^g and decay constants j^g take 
their values in the medium. 

The following combinations of probabilities are of interest to us, because if the time 
dependence could be measured and the different time behaviours disentangled, they would 
provide information on the effects that we are trying to uncover. To first order in perturbation 
theory we find 

P(K° -> K°) - P(K° -> K°) = 2Re5 K e _7xi - 2Re8 K e'^ st + Alm8 K sin(Amt) e -^ 7i . (21) 



[I]) does not depend on £33. It allows to measure both ReciK and Im5x CPT breaking 
parameters. 

P(K° -> K°) - P(K° -> K°) = -2Ree e _7L * - 2Ree e~^ st + 4Ree cos(Am i)e^ 7 ' .(22) 
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The coefficient of this combination of probabilities does not depend on the medium and 
would allow to extract Re e. 

The following sum singles out the parameters £33 and x and washes out the vacuum 
effects, 

P(K°^K )+P(K°^K ) = (I + x 33 )e-^ + (^-X33)e-^ 

+ (cos(Amt) + 

xsin(Amt))e~37*. (23) 

Thus we see that from observation of modulated oscillations we could in principle extract 
both the CPT violating and the decoherence parameters from the experiment. 

Finally, we consider the decays into two pion states from which e'/e has been recently 
measured through the double ratio 

T(K L -»■ TT + ir-)/T(K S -> 7T + 7T-) 



T(K L -> n°n )/T(K s -> 7r°7r°) 
In presence of matter, in the large time regime where only the long component survives, we 



obtain for it (compare with for the case of quantum gravity). 

1 + 6Re 



e' 







2 




2 +X 33 



Here we have used the definitions in (|J) replacing K s , K L and bu Ifg, and e^; A/g 
stand for the short and long components inside matter. Since X33 is positive, we see that the 
effect of entanglement with the environment would decrease the signal for e', i.e., it would 
be bigger than measured, e' has been measured in the vacuum whereas sl can be measured 
as well in the experiment. These decay modes thus provide also a way to determine X33. 



5. Conclusions 

We draw the attention to effects of interaction with the environment in neutral kaon envo- 
lution, that maybe of relevance in the analysis of the CP breaking parameters, made from 
more complete and precise data that will be available soon, in the future 0. The Lindblad 
master equation provides a universal form that incorporates the tiny corrections of non- 
hamiltonian evolution. On the other hand, decoherence represents a remarkable quantum 
phenomenon by itself. Thus we propose to share the experimental time for measurements of 
the decoherence parameter X33 in the medium, which we believe can be done by observation 
of time evolution. 
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7. Appendix 



In order to interpret the Lindblad piece of the evolution equation let us consider a simplified 
model of environment which consists only of two levels \E), \E'), with energies E' > E. 

The interaction kaon-environment is dominantly strong; it is dim, though, if the medium 
is diluted and the kaons are not too energetic. (We neglect the weak interaction of the 
kaons with the medium). The total hamiltonian has several terms : H = H + R + V. 
H Q = M I <g> I mv , R = I <g> (E\E)(E\ + E'\E')(E'\). V contains the interaction of the 
kaons with the environment; we do not include matrix elements that induce transitions like 
\Kq, E) — > \Kq, E), etc, diagonal in the kaons and in the medium. Having ordered the basis 
as \K ,E), \K ,E), \K ,E'), \K ,E'), V reads 

V = 




where A is a 2 x 2 matrix. The evolution equation 

ip(t) = [H,p(t)], 

reads 

ipi(t) = Wi(t),pi(t)], 

where 

pi (t) = e Hn +R)t p ty e -HH 0+ R)t 

V7(t) = e i(,H +R)ty e -i(H +R)t _ e iRty & -iRt _ 

The differential equation is equivalent to 



which leads to 



Pl {t) = Pl (0) - i f ' du\y!{u), Pl {u)\, 
Jo 

i Pl (t) = [W),P/(0)] -i fdulVrit),^)^^)]]. 

Jo 



Now we make the Markovian approximation, valid when the domain of integration that gives 
the dominant contribution is u — > t, so that the replacement pi(u) — > pi(t) is allowed. It 
amounts to neglecting memory effects. Furthermore we assume that the initial conditions 
factorize pi(0) = er(0) <E> p e nv, with a static environment p env = p\E)(E\ + p'\E')(E'\. 

We want to obtain an equation that takes into account the effects of the environment 
effectively, after having traced over the environment degrees of freedom: tr env pj(t) = cr(t). 

Given that if V = 0, Pi(t) = a{t) <E> p en v, w e expect pi(t) = a(t) <g> p env + 0(V). Also, 
tr env (Vj(t) Pl (0)) = tr env ( Pl (0)Vj(t)) = 0. 

Finally, we obtain 

ia{t) w Re f(t) [p'AA 1 + pA f A, <r(t)]+ 
+ilm f(t) (p(2AaA t - a A* A - A^Aa) + p '(2AVA - aAA* - AAV)) , 

where 

-i{E'-E)t _ I 

m = E>-E • 
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and we have only kept terms up to 0(V 2 ) included. The first term gives a correction 
to hamiltonian evolution and the second is a Lindblad like piece. Since our model contains 
only two levels for the environment these contributions depend on time, through the function 
f(t). It is only in the situation that the environment has a continuum spectrum of energies 
that the limit of large t gives, in the sense of distributions 

Im f(t) — > tv5(E' — E); 

it only contributes when there is a nonzero density of environment states, degenerate with 
the initial state. 

The Lindblad operators are the matrix elements A and their conjugates that, under 
evolution, would produce entanglement of the kaons with the environment. It also depends 
on the energy distribution in the environment. 
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